Inelastic current-voltage characteristics of atomic and molecular junctions 
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We report first-principles calculations of the inelastic current- volt age (I-V) characteristics of a 
gold point contact and a molecular junction in the nonresonant regime. Discontinuities in the I-V 
curves appear in correspondence to the normal modes of the structures. Due to the quasi-one- 
dimensional nature of these systems, specific modes with large longitudinal component dominate 
the inelastic I-V curves. In the case of the gold point contact, our results are in good agreement with 
recent experimental data. For the molecular junction, we find that the inelastic I-V curves are quite 
sensitive to the structure of the contact between the molecule and the electrodes thus providing a 
powerful tool to extract the bonding geometry in molecular wires. 



Inelastic scattering between electrons and phonons in 
a current-carrying wire is a source of energy dissipation 
for electrons. However, it can also yield a lot of informa- 
tion on the underlying atomic structure of the wire. This 
information can be extracted indirectly from the disconti- 
nuities in conductance that occur when the external bias 
is large enough to excite discrete vibrational modes of the 
wire. 0[ Recent experimentsjDn transport properties of 
atomic 2j and molecular 3, 4] junctions have indeed re- 
vealed such inelastic features. It is, however, not straight- 
forward to relate these features to specific vibrational 
modes. A nanoscale junction (often described as a quasi- 
one-dimensional system) with N atoms supports 3N vi- 
brational modes. In a strictly one-dimensional system, 
only longitudinal modes can be excited via electronic 
coupling. However, the modes of a realistic junction are 
not necessarily purely transverse or purely longitudinal 
with respect to the direction of current flow, j^, 0, 
Therefore, the inelastic current- volt age (I-V) character- 
istics are likely to depend strongly on the detailed atomic 
structure of the full system. This is particularly relevant 
for molecular junctions for which the contact geometry 
between the molecule and the bulk electrodes is difficult 
to control in experiments. |^ |^ 

In this letter we first derive an expression for the inelas- 
tic current in a current-carrying system in terms of scat- 
tering wavefunctions. This expression allows us to study 
the inelastic I-V characteristics of a given nanoscale junc- 
tion using first-principles approaches. As an example we 
study the effect of vibrations on the electron dynamics in 
a gold point contact and a single- molecule junction. For 
the gold point contact, the magnitude of the calculated 
inelastic current as well as its onset compare very well 
with recent experimental results. 2] For the molecular 



junction, we analyze the case in which the molecule is 
equally bonded to the two bulk electrodes and the case 
in which the two contacts are different. We find that 
the inelastic I-V characteristics are very different in the 
two cases. This result shows that inelastic spectroscopy 
could be used quite effectively to extract information on 
the contact geometry of molecular wires. 

Let us start by deriving an expression for the inelas- 
tic current. We assume that the phonon distribution is 
at equilibrium at all (small) biases, thus neglecting lo- 
cal heating. 5] We have previously shown that for small 
biases this effect is small, provided good thermal con- 
tacts exist between the nanostructures and the bulk elec- 
trodes. The many-body Hamiltonian of the system is 
(atomic units are used throughout this paper) [5| 
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where Hei is the electronic part of the Hamiltonian; 
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bution where is the normal coordinate and uOi^ is the 
normal frequency corresponding to the i-ih. ion and ji- 
th component; finally, Hei-vib describes the electron- ion 
interaction and has the following form: 
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where a = L^R; and bjjy are the electron and phonon 
annihilation operators, respectively, satisfying the usual 
commutation relations. Ai^jj^ are the matrix elements of 
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the transformation from cartesian coordinates to normal 
coordinates, and J^EuE2 electron-phonon coupling 

constant which can he directly calculated from the scat- 
tering wave-functions 
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where we have chosen to describe the electron-ion in- 
teraction with pseudopotentials V^^ (r, Ri) for each i-th 
ion. 



Similar to what has been done in Ref. 6, we treat 
the electron-phonon interaction to first-order perturba- 
tion theory. pj| Due to the orthogonality condition be- 
tween phonon states, higher harmonics for each phonon 
mode appear only in third-order perturbation theory 
and are therefore small. 12] We develop the full many- 



body wavefunctions in terms of the states ( ^ 
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{njiy\. The single-particle electronic state is 

described by ^^^^^ (r, Ky), corresponding to electrons 
incident from the left (right) electrodes with energ y E 
and momentum Ky parallel to the electrode surface, [l^ 
These electronic states are calculated self-consistently by 
means of a scattering approach within the density func- 
tional theory of many-electron systems. 10] The phonon 
state is described by {njjy\, where rijjy is the number of 
phonons in the ju-th normal mode. 
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where (n^j,) = 1/ [e-x.-^ {fujOjj, /kBTy^) — 1] is the Bose- 
Einstein distribution per mode at a given wire temper- 
ature Tu;, and () indicates the statistical average. The 
above expression allows us to calculate the inelastic cur- 
rent. It is evident from Eq. O that for a fixed partial 
density of states, the magnitude of the inelastic current 
is determined by the coupling constant J^Ei^E2 
transformation matrix A = {Ai^jj^} which contains the 
information on the geometry of the structure and hence 
on the character (transverse versus longitudinal) of the 
different modes. 14] We will be concerned with the ex- 
tra inelastic current due to the vibrational modes of the 
atoms of the nanoscale constriction with respect to the 
continuum spectrum of modes of the bulk electrodes. 
If the electronic temperature Tg is zero, then, for an ex- 
ternal bias V, only those normal modes with eigenener- 
gies hujjiy < eV can be excited and contribute to Eq.[5J In 
addition, due to our assumption of negligible local heat- 
ing, the averaged number of phonons (n^v) is zero for all 
normal modes. In this case the first-order correction to 
the current induced by electron-phonon interaction as- 
sumes the following simple form: 
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In the above expression is the partial density of 
states corresponding to the current-carrying states 
and e{E^njiy) = E -\- {rijy + 1/2) hujjiy is the energy of 
state 1^^; n^v). We have also assumed that the electrons 
rapidly thermalize into the bulk electrodes so that their 
statistics are given by the equilibrium Fermi-Dirac dis- 
tribution, f^^"^^ = l/(exp[(£; - EFLiR))/kBTe] + 1) with 
chemical potential £^FL(i?) deep into the left (right) elec- 
trode. [13] Using lime^o = ^(l) + i7rS{z), the first- 
order correction |(5^^;njj,) assumes the following form: 
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where only the left-travelling electronic states contribute 
(if the left electrode is positively biased). 

We are now ready to use the above expression to study 
inelastic scattering in specific systems. We choose to first 
study a gold point contact for which experimental results 
are available ^% and then discuss the case of a molecular 
junction. In Fig. Q we plot the inelastic conductance for 
a single gold atom. 

In the absence of inelastic scattering and for the bias 
range of Fig. ^ the I-V characteristics of this system are 
linear with differential conductance G 1.1 Go, where 
Go = 2e^ /h. 5] When electron-phonon interactions are 
considered, two transverse modes with energy hw :^ 10.8 
meV are first excited with increasing bias. However, 
due to their transverse character, these modes contribute 
negligibly to the inelastic current as determined by the 
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FIG. 1: Absolute value of the differential conductance due to 
electron-phonon interaction as a function of bias for a gold 
point contact. Two normal modes corresponding to trans- 
verse vibration of the gold atom between the electrodes have 
energies of 10.8 meV, which are close in energy to the longi- 
tudinal mode at 11.5 meV (shown in the figure). 



product between the transformation matrix A and the 
coupling constant J^'^. An abrupt change in differen- 
tial conductance appears at V 11.5 mV corresponding 
to the excitation of a longitudinal vibrational mode (see 
schematic in Fig. 16] Both the onset bias as well as 
the change in conductance (of about 1%) are in good 
agreement with experimental reports on gold point con- 
tacts. The longitudinal and transverse modes are very 
close in energy but mainly the longitudinal mode con- 
tributes to the inelastic current, j^, [2] so that in exper- 
iments the transverse ones would not be easily resolved. 
This is even more evident in the case of the molecular 
junction. 

In Fig. [21 we plot the inelastic conductance in the case 
in which a phenyldithiolate molecule forms symmetric 
contacts on both sides of the junction, i.e., each S atom 
is bonded to a fiat surface. In this case there is a total 
of 14 modes with energy less than 100 meV. A promi- 
nent change in conductance occurs at a bias of about 18 
mV, i.e., at a bias large enough to excite two modes with 
large longitudinal component (see Fig. ^3 The in- 
elastic contribution from two transverse modes at lower 
bias jSi] is almost four orders of magnitude smaller. Simi- 
larly, three quasi-transverse modes with energies between 
20 and 50 meV contribute negligibly to the inelastic con- 
ductance. They only appear as small features in the sec- 
ond derivative of the current with respect to the bias (see 
Fig. [2|). 18] It is likely that due to noise and other ef- 
fects, such modes would not be resolved in experiments. 
Increasing the bias further, a second large step in the 
absolute value of the conductance is found at about 50 
mV (see Fig. 0. This again corresponds to a predomi- 
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FIG. 2: Absolute value of the differential conductance due to 
electron-phonon interaction as a function of bias for a sym- 
metric molecular junction. The derivative of the conductance 
with respect to bias is also shown (a broadening of 1 meV 
is introduced). The schematics show only the modes that 
contribute the most to the inelastic current. 



nantly longitudinal mode (see schematic in Fig.I2|). This 
mode is then followed by others that have both a trans- 
verse and a longitudinal component. The magnitude of 
the conductance steps depends on the relative amount of 
the two components as well as the product between the 
transformation matrix A and the coupling constant. ^3 

We conclude by showing how sensitive the inelastic cur- 
rent is to any change in the bonding properties of the 
molecule to the electrodes. We illustrate this in Fig. IHl 
where we plot the inelastic conductance for the same 
molecule but with one of the S atoms bonded to a H 
atom which, in turn, is not bonded to the nearby surface 
(see schematics in Fig.|2|. 19] Such a configuration could 
be easily realized in experiments. Isl In the present case 
there are 13 modes below 100 meV, of which only six 
have large longitudinal component (shown in Fig. 
and contribute to small steps in the inelastic conduc- 
tance, with the one at about 11 meV showing the largest 
relative contribution. All other modes contribute neg- 
ligibly to the current. Comparing Figs. [21 and 01 it is 
clear that contact geometry affects considerably the I-V 
characteristics of the system both in the position of the 
inelastic discontinuities as well as in their relative mag- 
nitude. This fact, however, can be used advantageously 
to extract a posteriori the contact structure of molecular 
junctions thus providing a powerful diagnostic tool for 
nanoscale electronics. 
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FIG. 3: Absolute value of the differential conductance due to 
electron-phonon interaction as a function of bias for a molec- 
ular junction with asymmetric contacts. The derivative of the 
conductance with respect to bias is also shown (a broadening 
of 1 meV is introduced). The schematics show only the modes 
that contribute the most to the inelastic current. 

Biomedical Institute. Acknowledgment is also made to 
the Donors of The Petroleum Research Fund, adminis- 
tered by the American Chemical Society, for partial sup- 
port of this research. One of us (MZ) acknowledges par- 
tial support from an NSF Graduate Fellowship. 



E-mail address: diventra@physics.ucsd.edu 
^ N. Lorente, M. Persson, L.J. Lauhon, W. Ho, Phys. Rev. 

Lett. 86, 2593 (2001). 
^ N. Agrait, C. Untiedt, G. Rubio-Bollinger, and S. Vieira, 

Phys. Rev. Lett. 88, 216803 (2002). 
2 H. Park et al. Nature, 407,57 (2000); L.H. Yu and D. Na- 

telson, Nano Letters, 4, 79 (2004). 

R.H.M. Smit, Y. Noat, C. Untiedt, N.D. Lang, M.C. van 
Hemert, and J.M. van Ruitenbeek, Nature 419, 906 (2002). 
^ Y.-C. Chen, M. Zwolak, and M. Di Ventra, Nano Lett. 3, 
1691 (2003). 

^ M.J. Montgomery, J. Hoekstra, T.N. Todorov, and A. Sut- 
ton, J. Phys.: Cond. Mat. 15, 731 (2003). 
A. Troisi, M.A. Ratner, and A. Nitzan, J. Chem. Phys. 
118, 6072 (2003). 

^ X. Xiao, B. Xu, and N.J. Tao, Nano Lett. 4, 267 (2004). 

^ M. Di Ventra, S.T. Pantelides, and N.D. Lang, Phys. Rev. 
Lett. 84, 979 (2000). 



^° N.D. Lang, Phys. Rev. B 52, 5335 (1995); M. Di Ventra 
and N.D. Lang, Phys. Rev. B 65, 045402 (2002); Z. Yang, 
A. Tackett, and M. Di Ventra, Phys. Rev B. 66, 041405 

(2002) . 

Since we are in the regime of weak electron-phonon cou- 
pling we neglect the elastic phonon renormalization of the 
electronic spectrum (see e.g. G.D. Mahan, Many Particle 
Physics, 2nd Edition (Plenum Publishers, NY, 1990), p. 
285). 

^2 V. Aji, J.E. Moore, and CM. Varma, |cond-mat/0302222( 
A. Mitra, I. Aleiner, A.J. Millis, cond-mat/0311503 
Y.-C. Chen and M. Di Ventra, Phys. Rev. B 67,153304 

(2003) . 

A quantitative way of distinguishing between longitudinal 
versus transverse character of each mode is to sum up the 
absolute values of all the z-component (i.e., along the cur- 
rent flow) matrix elements of A. 

We have employed Hartree-Fock total energy calculations 
[see, e.g., J. A. Boatz and M. S. Gordon, J. Phys. Chem., 
93, 1819 (1989)] to evaluate the vibrational modes of the 
single- molecule junctions and the gold point contact. For 
these calculations, the assumed surface orientation is [111] 
represented by a triangular pad of gold atoms with infi- 
nite mass (see insets of Figs. Q] and 0). The structures were 
relaxed starting from an initial geometry where the single 
gold atom of the point contact faces the center of the trian- 
gular pad at a distance of 2.3 A; and for the single- molecule 
junction with symmetric contacts the initial S-surface dis- 
tance is 2.4 A. In the case of asymmetric contacts the S 
atom on one side of the junction is bonded to a H atom 
(see insets of Fig. This thiol termination is assumed 
to be far from the nearby surface so that the modes of 
this structure are not influenced by the gold atoms of the 
surface. 

Note that the inelastic correction to the current is actually 
negative. For convenience we plot in this paper the absolute 
value of this correction. 

In the molecular junction cases, the dominant net contri- 
bution to the electron-phonon coupling is from the S atoms 
and the adjacent C atoms (i.e., the two C atoms that form 
a straight line with the S atoms). This implies that those 
modes that do not have large displacements of these atoms 
along the z direction will contribute negligibly to the in- 
elastic current, even if the global character of the mode 
is longitudinal as inferred from the procedure outlined in 
Ref. [13 . This is, e.g. the case for the mode at about 90 
meV shown schematically in Fig. |21 

A broadening of 1 meV is introduced in all calculations 
to make the second derivative of the inelastic current as a 
function of bias finite at the conduction jumps. 
In order to facilitate the comparison with the previous case, 
we have assumed in this calculation the same coupling con- 
stants, partial density of states and unperturbed current 
as for the symmetric junction case. 



